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Abstract:  We  report  observation  of  query  pulse  length  dependent  Ramsey 
interference  (QPLD-RI),  using  pulsed  Raman  excitation  in  rubidium  vapor. 

This  is  observed  when  a  long,  attenuated  query  pulse  is  used  during  pulsed 
Raman  excitation.  We  explain  the  physical  mechanism  behind  the  QPLD- 
RI  using  a  Bloch  vector  model.  We  also  use  numerical  solutions  to  time- 
dependent  density  matrix  equations  to  simulate  this  interference  effect, 
showing  qualitative  agreement  with  experimental  results.  Presence  of  such 
interference  could  create  a  potential  source  of  error  in  a  vapor  cell  Raman 
clock  constructed  using  frequency-domain  Ramsey  interference  (FDRI). 
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1.  Introduction 

Optical  (or  microwave)  fields  separated  in  space  (or  time)  can  interfere  to  produce  Ramsey 
fringes  via  the  coherence  induced  by  the  atomic  polarization  [1].  During  such  a  process,  the 
coherence  is  transported  by  an  atom  (or  molecule)  and  interference  is  produced  by  the  read¬ 
out  field  depending  on  the  phase  difference  between  the  oscillating  field  and  the  coherence. 
Over  the  years,  several  papers  have  reported  Ramsey  interference  using  different  schemes 
such  as  the  Doppler-free  two-photon  resonance,  ground-state  Zeeman  coherence  and 
polarization  rotation,  and  suggested  their  applications  in  precision  metrology  [2-6].  Some  of 
these  studies  have  also  shown  interference  resulting  from  atomic  diffusion  where  a  coherently 
prepared  atom  escapes  and  reenters  the  interaction  region  with  a  residual  phase  memory  [5,6]. 

Frequency-domain  Ramsey  fringes  are  generated  by  changing  the  frequency  of  the 
oscillatory  field  away  from  the  atomic  resonance  during  interrogation.  This  method  was 
originally  applied  to  separated  field  microwave  interactions  in  an  atomic  beam  to  yield 
narrower  atomic  resonances  [7].  It  is  now  widely  used  in  the  development  of  cesium  beam 
atomic  clocks  [8,9].  Later,  separated-field  Raman  interaction  using  a  three-level  A-system 
was  developed  to  replace  the  microwave  interaction  in  the  atomic  beam  by  an  all-optical 
interaction  while  generating  Ramsey  interference  [10-12].  The  interference  produced  by  this 
scheme  is  commonly  referred  to  as  Raman  Ramsey  interference  (RRI). 

Recently,  an  extension  of  this  approach  to  time-domain  has  been  studied  in  atomic  vapor 
cells  with  renewed  interest  [13-16].  This  technique  can  be  applied  to  realize  a  potentially 
high-performance  vapor  cell  Raman  clock  compared  to  that  using  coherent  population 
trapping  (CPT)  [17,18].  In  a  recent  study,  we  have  shown  that  high-contrast  frequency- 
domain  Ramsey  interference  (FDRI)  fringes  with  sub-kilohertz  linewidth  can  be  produced  in 
coherently  prepared  85Rb  atomic  vapor  [19].  A  repeated  sequence  of  time-separated  optical 
Raman  pulses  was  used  to  generate  FDRI  fringes. 

In  this  paper,  we  show  a  new  type  of  Ramsey  interference  effect  using  a  similar  pulsed 
excitation  mechanism  in  rubidium  vapor.  Unlike  FDRI,  these  fringes  are  observed  as  time- 
domain  oscillations  in  the  transmission  amplitude  of  a  long,  attenuated  query  pulse.  As 
discussed  later,  the  frequency  of  these  oscillations  is  inversely  related  to  the  Raman  detuning, 
and  the  number  of  observed  oscillations  depends  on  the  duration  (or  length)  of  the  query 
pulse.  We  refer  to  them  as  query  pulse  length  dependent  Ramsey  interference  (QPLD-RI). 
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We  use  a  Bloch  vector  model  to  provide  a  clear  physical  explanation  of  RRI  observed  in 
the  form  of  FDRI  and  QPLD-RI.  We  also  use  a  numerical  model  to  describe  the  dynamics  of 
the  A-system  under  pulsed  Raman  excitation,  agreeing  closely  with  our  experimental  results. 
Keeping  in  mind  a  possible  high-performance  vapor  cell  clock  development,  we  have  also 
explored  the  magnetic  field  dependence  of  FDRI  and  QPLD-RI  using  a  clock  transition  in 
85Rb  atom.  RRI  can  also  be  used  in  sensitive  magnetometry  [20]  if  Zeeman  sublevels  with 
first  order  magnetic  field  sensitivity  are  used  in  Raman  transitions  during  pulsed  Raman 
excitation. 


RP-1  RP-2 
(C01/C02)  (g>i/G>2) 

i  \ 

T 


(b) 

Fig.  1.  (a)  Raman  excitation  in  a  three-level  A-system:  coi  and  co2  correspond  the  frequencies  of 
the  Raman  beams,  A  and  S  correspond  to  the  difference  (or  two-photon)  and  average  frequency 
detuning,  respectively  (b)  repeated  sequence  of  Raman  pulse  (RP)  used  in  generating  Ramsey 
interference. 

Our  study  shows  that  QPLD-RI  can  produce  pulse  length  dependent  broadening  as  well  as 
asymmetry  in  FDRI  fringes  used  for  a  vapor  cell  Raman  clock,  causing  a  potential  source  of 
error.  Specifically,  we  show  that  broadening  affects  the  fractional  frequency  stability  of  a 
clock,  and  asymmetry  changes  the  magnitude  of  the  light-shift  induced  error  due  to  optical 
detuning,  which  is  unavoidable  in  a  vapor  cell  based  Raman  clock  because  of  the  presence 
several  hyperfine  transitions  in  the  excited  state  manifold,  which  all  contribute  to  the  clock 
signal  because  of  Doppler  broadening.  Thus,  in  choosing  the  optimal  parameters  for  such  a 
clock,  the  effect  of  QPLD-RI  has  to  be  taken  into  account. 

2.  Theoretical  description 

Two-photon  resonant  Raman  excitation  in  a  three-level  A-system  forms  the  basis  for  both 
QPLD-RI  and  FDRI  to  be  discussed  here.  Figure  la  shows  such  a  scheme.  Here,  &>7  and  co2 
are  the  frequencies  of  the  two  excitation  fields,  detuned  from  resonance  by  Sj  and  S2 , 
respctively;  £?7  and  Q2  correspond  to  the  Rabi  frequencies  for  the  two  legs;  fe(<57+<52)/2  is  the 
average  detuning  of  the  two  fields,  and  A ={SrS2)  is  the  difference  detuning,  which  is  also 
called  the  two-photon  detuning  or  the  Raman  detuning.  We  assume  that  the  excited  state 
population  decays  equally  to  the  two  metastable  ground  states,  with  a  net  decay  rate  of  T.  If 
the  Raman  beams  are  continuous,  CPT  resonance  can  be  observed  in  the  transmitted  laser 
field  [21-24].  On  the  other  hand,  if  the  atomic  medium  is  excited  by  a  repeated  sequence  of 
pulsed  Raman  beams  as  shown  in  Fig.  lb,  RRI  can  be  observed.  In  this  case,  a  relatively  long 
pulse  (duration  Xi)  is  used  first  to  prepare  the  atoms  in  the  CPT  state,  also  known  as  the  dark 

state  (|-)  =  (Q2  |l)-Q1  |2))/^/qi2  +Q2  )•  After  a  free  evolution  time,  T,  another  Raman  pulse 

of  short  duration  x2  is  used  as  a  detection  pulse  to  observe  Ramsey  interference.  The 
interference  is  created  by  the  second  pulse  via  the  atomic  coherence  induced  by  the  first 
Raman/CPT  pulse.  Short  detection  (or  query)  pulse  is  generally  used  in  order  to  avoid  any 
significant  perturbation  of  coherence  produced  by  the  first  pulse  [14].  Choosing  a  longer  x2 
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will  tend  to  suppress  the  interference  effect  unless  one  uses  a  strongly  attenuated  query  pulse, 
as  discussed  in  greater  detail  later.  Experimentally,  the  interference  fringes  are  observed  by 
measuring  the  transmitted  amplitude  of  the  short  detection  pulse  using  electronic  gating  and 
time  integration  while  changing  the  interference  condition  by  changing  the  difference 
frequency,  A,  between  the  Raman  beams.  To  gain  further  insight  into  the  origin  of  Raman 
Ramsey  interference  phenomena  and  their  characteristics,  here  we  will  use  a  three-element 
Bloch  vector  (BV)  model  [22]  to  describe  the  time  evolution  of  the  A-system  under  pulsed 
Raman  excitation  and  explain  the  formation  of  FDRI  and  QPLD-RI. 

If  we  consider  the  decay  rate  T  of  the  excited  state  to  be  much  greater  than  any  other  rates 
in  the  A-system,  the  unpopulated  atomic  state  |2)  can  be  eliminated  from  the  interaction 

picture  using  an  adiabatic  approximation.  Such  an  approximation  is  fairly  consistent  with  our 
experimental  study  of  interest.  In  this  case,  the  time  evolution  of  the  A-system  can  be 
represented  by  the  motion  of  a  three-element  B  V.  The  components  of  B  V  can  be  described  by 
the  density  matrix  elements  of  the  atomic  states  |l)  and  |3)  alone.  The  BV  in  this  case,  is 

somewhat  analogous  to  the  BV  of  a  two-level  system,  and  has  been  discussed  at  length  for 
Raman  excitation  in  ref.  22. 


A  partially  diagonalized  basis  consisting  of  the  dressed  states  +) 


1 

12 


(|1)  +  |3))  for  Q, 


-  Q2  =  Q  provides  the  simplest  interpretation  for  the  BV  model.  The  equation  of  motion  for 
the  Raman  BV  vector  [22]  is  given  by 


R-QxR  -  a  R  +  a  e  3 


(1) 


where  R  =  RX  ex  +  R2e2  +  R3  e3  is  the  BV,  and  Q  =-  Ae1  +  p  e3  is  the  so-called  torque  vector. 
The  three  real  elements  (Rlf  R2,  R3)  of  the  BV  are  defined  in  terms  of  the  density  matrix 

elements  in  the  dressed  basis  as  p_+  =  ^(Rx+iR2) ,  and  p__-  p++  =  R3 .  The  parameter 

a  -  Q^/2T  =  Q2/2T  corresponds  to  the  Raman  pumping  rate,  fi  =  fl2R  S/ (r2  +4S2^J 
corresponds  to  the  frequency  difference  between  |  +)  and  |  -)  states  which  in  our  discussion 

is  generally  considered  to  be  negligibly  small  (i.e.  /?  -  0 ),  and  C1R  +  £l22  is  the 

Raman  Rabi  frequency.  Thus,  the  motion  of  the  BV  is  primarily  governed  by  the  choice  of 
parameters  a  and  A.  For  a  closed  atomic  system,  p__  +  p++  =1 .  Note  that  the  Raman  BV  does 


not  conserve  its  amplitude  at  all  times  due  to  the  presence  of  decay  of  R  and  growth  of  R3 
(second  and  third  terms  in  Eq.  (1).  As  can  be  seen  from  Eq.  (1),  it  can  start  from  zero 
amplitude,  and  grow  into  a  maximum  steady-state  amplitude  of  one.  This  is  to  be  contrasted 
with  the  case  of  a  two  level  system  without  decay  (e.g.,  a  microwave  transition  under  ideal 
conditions).  In  such  a  system,  the  length  of  the  BV  is  conserved.  However,  if  decay  is 
included  in  a  two  level  system  (e.g.,  an  optical  transition  where  the  upper  level  decays 
rapidly),  the  length  of  the  BV  is  also  not  conserved,  even  if  the  system  is  closed. 

Note  that  the  BV  described  above  can  also  be  expressed  in  the  atomic  basis  as 


R  =  R 1  ex  +  R  2  e2  +  R  3  e2 ,  where  pX3  =  —  ^Rx  +  i  R  2 )  and  pn  -  p33  =  R3 .  The  ei  basis  is  related 


to  the  ei  basis  by  a  rotation  of  angle  26  in  the  1-3  plane,  where  tan  6  -  Qj  /  Q2  [22].  For  the 


case  of  f27  =  Q2 ,  the  angle  of  rotation  is  26  =  n  /  2  ,  so  that  Rl  =  R3,  R2  =  R2,  R3  =  -R{  .  This 
corresponds  to  the  following  explicit  relations  between  the  elements  of  the  density  matrix  in 
these  two  bases:  Re (pl3 )  =  - ( p__  -  p++ )  /  2  ,  hn  (p13 )  =  Im ( p_+ ) ,  and  pn-  p33  =  2 Re ( p_+ ) . 
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However,  the  dynamics  of  the  system  is  more  transparent  in  the  dressed  states  basis  [22]; 
hence,  we  use  this  basis  throughout  this  paper. 


ft 


Fig.  2.  Bloch  vector  model  showing  the  dynamics  of  a  three-level  A-system  during  pulsed 
Raman  excitation.  The  figures  show  the  BVs  (a)  after  interaction  with  the  first  Raman/CPT 
pulse  (dotted  blue  arrow:  start  of  BV;  solid  blue  arrow:  end  state  of  BV)  (b)  during  free 
evolution  time,  T,  and  (c)  after  interaction  with  the  second  Raman  (or  detection)  pulse. 

The  rotation  of  the  Raman  BV  during  interaction  with  the  pulsed  Raman  beams,  and 
during  free-time  evolution,  T,  can  be  described  as  follows.  The  first  CPT  pulse  (shown  in  Fig. 
lb)  causes  a  strong  interaction  due  to  higher  pulse  intensity  and  long  ii.  Under  this  condition, 
we  can  assume  that  a  »  |A|  and  arx  »  1 .  This  causes  the  system  to  reach  the  steady  state 
after  interaction  with  the  first  pulse.  The  steady  state  solution  to  the  equation  of  motion  in  Eq. 
(1)  is  found  from  R  =  0 ,  which  indicates  formation  of  a  complete  dark  state  described  by  the 
BV  l.e3  with  R3  =  1 .  This  also  corresponds  to  p__  -  1,  p++  =  0,  p_+  =  0  .  Figure  2a 

shows  pictorially  evolution  of  the  BV  starting  from  R(0)  ,  indicated  by  the  dotted  blue 
arrow,  to  /^(rj  ,  indicated  by  the  solid  blue  arrow,  during  interaction  with  the  first  pulse,  for 
\P\  » |A|  .  Under  this  condition,  the  torque  vector  is  Q  «  fie3 ,  which  causes  the  spiraling 
motion  around  the  e3  axis.  For  the  case  of  P  =  0 ,  the  BV  simply  moves  to  the  steady  state 
value  in  a  straight  line,  starting  from  any  initial  value  [22].  To  illustrate  the  motion  of  the  BV 
in  the  free  evolution  zone  and  the  query  pulse  zone,  we  assume  p  -  0  in  these  zones. 

In  the  absence  of  interaction  during  the  free  evolution  time  T,  a  =  0.  Under  this  condition, 
one  can  no  longer  ignore  the  two-photon  detuning,  A,  which  causes  the  Bloch  vector  R  to 
precess  at  a  rate  A  around  the  ex  -axis,  as  shown  in  Fig.  2b.  Using  Eq.  (1),  one  can  find  the 

Raman  BV  at  time  t  =  rx+T  as  R(tx+T}  =  sin(A .T)e2  +cos(A T)e3  .  This  implies  that  at 
the  end  of  time  T,  the  BV  sweeps  a  total  angle  0  (  =  A.T).  For  example,  if  A  =  0 ,  the  BV  does 
not  rotate  and  R[rx+T^-\.e3  .  During  interaction  of  the  atoms  with  a  repeated  sequence  of 
Raman  pulses,  A  can  be  changed  slowly  around  A  =  0,  thus,  causing  the  BV  R{tx+T^  to 
oscillate  in  the  [e2  -e3)  plane.  The  effect  of  this  oscillation  can  be  observed  by  measuring  the 
absorption  of  a  second  Raman  (or  detection)  pulse  which  is  present  during  the  time 
t  =  tx+T  and  t  =  rx  +  T  +  t2  . 
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The  BV  at  t  =  rl+T  serves  as  the  initial  BV  for  interaction  with  the  second  Raman  pulse. 
It  is  assumed  that  the  optical  fields  in  the  second  Raman  pulse  have  the  same  phase  as  that  of 
the  first  Raman  pulse  so  that  the  motion  of  the  B  V  simply  continues  from  where  it  was  left  off 
at  the  end  of  time  t  =  zx  +  T  .  During  this  interaction  in  the  query  zone,  we  also  assume 

a  »  |A|  in  Eq.  (1),  and  the  motion  of  the  BV  is  described  as 
R{rx  +T  +  t2)=  sin(AT)  e~aTl  e2  +  |^l-(l-cos(AT)  e~aTl  e3  .  While  observing  FDRI,  z2 

is  kept  short  so  that  a  r2 « 0 ,  and  the  Raman  B  V  practically  remains  unchanged  after 
interaction  with  the  second  pulse.  This  effect  is  shown  in  Fig.  2c.  It  is  obvious  that  if  x2  is 
made  long,  the  BV  will  eventually  grow  linearly  to  the  steady  state  BV  given  by  Rs  -\.e3  , 
indicating  the  formation  of  a  new  dark  state.  The  transmitted  amplitude  of  the  second  pulse 
directly  corresponds  to  the  projected  BV  component  R3(  -  cos (AT)  ),  which  is  also  a 

measure  of  the  dressed  state  population  difference  (/?_  _  -  p++ ) .  FDRI  is  observed  as 

oscillations  in  R3  caused  due  to  change  in  A  during  repeated  interaction  with  the  pulse 

sequence.  The  frequency  width  of  these  fringes  is  inversely  related  to  T,  and  is  not  affected  by 
power  broadening. 


Fig.  3.  Plots  showing  oscillations  in  BV  components  R2  and  R3  as  a  function  of  x2  for  A  =  10 
KHz,  a  =  0.1  A  and  T  =  100  ps. 

The  BV  model  also  explains  the  new  type  of  Ramsey  interference,  QPFD-RI.  This 
interference  is  observed  when  the  amplitude  of  the  second  Raman  pulse  is  significantly 
attenuated  and  the  duration  r2  of  the  pulse  is  made  long.  In  this  case,  the  motion  of  the  BV  in 
the  query  zone  is  governed  by  both  a  and  A.  Using  appropriate  initial  conditions  at 
t  =  rx  +T  ,  the  solution  to  Eq.  (1)  is  obtained  as 


*K)  = 


|^sin(A.(r2  +T7))  -  cos  6.  sin(A.(r2  +0))]  e 


sin  2  <9 


(2) 


|^cos(A.(r2  +T))  -  cos  6.  cos(A.(r2  +0))]  e 


t2  +- 


1  +  cos20 
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where  6  =  arccos^a/ ^a2  +  A2  j .  The  projected  BV  component  R 3  describes  the  Ramsey 

interference  observed  in  the  transmitted  amplitude  of  the  query  pulse  as  a  function  of  r2.  It 
contains  two  damped  oscillatory  terms  decaying  at  a  rate  of  a.  Figure  3  shows  the  oscillations 
in  BV  components  R2  and  R3  as  a  function  of  r2.  We  chose  A  =  10  KHz,  and  a  =  0.1  A  for 

this  case.  Also,  since  a  |A| ,  0  »90°,  in  this  case,  the  BV  can  be  approximated  as 

R(  r2)  -  sin(A.(r2  +  T))  e~aZl  e2  +  cos(A.(r2  +  T))  e~aTl  e3  .  One  can  draw  the  following 

conclusions  about  QPLD-RI.  Temporal  width  of  these  fringes  is  inversely  related  to  A.  The 
fringes  will  disappear  if  A  is  matched  to  the  two-photon  resonance  condition  (i.e.  A  =  0).  The 
initial  phase  of  the  oscillation  is  determined  by  T.  Decay  time  of  the  fringes  is  inversely 
related  to  the  Raman  pumping  rate  a  of  the  query  pulse.  The  present  BV  model  does  not 
include  ground-state  decoherence,  y  and  its  effect  on  the  decay  of  QPLD-RI. 

Even  though  this  model  provides  a  simple  understanding  of  the  dynamics  of  the  Ramsey 
interference,  a  quantitative  modeling  using  time-dependent  density  matrix  equations  is 
necessary  for  describing  these  interference  phenomena  more  accurately.  We  will  briefly 
describe  the  basic  density  matrix  equations  and  present  our  numerical  simulation  results  using 
such  a  model.  The  density  matrix  equations  for  a  three-level  atomic  system  are  described  by 
the  Liouville  equation: 


P  =  ~:j[H>p]  +  P**,  (3) 

where  p  represents  the  density  operator  in  the  rotating-wave  basis  consisting  of  nine  matrix 
elements  for  a  three-level  system.  pdecay  describes  various  possible  relaxation  terms  between 
the  atomic  levels.  The  diagonal  elements  of  p  correspond  to  the  population  terms  satisfying 
the  condition  pn  +  p22  +  p33  - 1  for  a  closed  system,  and  the  off-diagonal  elements 
correspond  to  the  coherence  between  various  atomic  levels  created  due  to  interaction  with  the 
laser  fields.  The  interaction  Hamiltonian  for  such  a  three-level  system  in  the  time  independent 
basis  is  given  by 
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where  the  parameters  are  as  defined  earlier  in  Fig.  1.  Interaction  of  the  atoms  with  the  Raman 
pulse  sequence  and  generation  of  Ramsey  interference  can  be  described  by  an  exact  solution 
to  eight  independent  time-dependent  equations  given  in  Eq.  (3).  We  obtained  the  numerical 
solutions  to  these  equations  during  three  sequential  time  intervals  xh  T,  and  r2  respectively, 
and  used  the  solution  from  one  time  interval  as  the  initial  values  for  the  next  time  interval. 

Figure  4a  shows  a  typical  example  of  numerically  computed  FDRI  fringes.  For  this  case, 
we  used  a  Raman  pulse  sequence  with  parameters  r7  =  4  x  103  -(l/T)  [»212  ps  for  r  =  2k  (6 
MHz)],  T  =  795  jus ,  and  r2  =  1  jus.  The  Rabi  frequencies  for  the  Raman  beams  were  taken  to  be 
equal,  and  weak  compared  to  the  excited  state  decay  rate:  Ql  -  Cl2- 27t(T / 10) ,  and  3  was  set 
equal  to  zero.  Amplitude  of  FDRI,  in  this  case,  was  calculated  by  finding  the  integrated 
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rx+T  +  r2 

signal  J  Im[/?13]-  dt  as  a  function  A.  The  line  width,  Av  of  the  Ramsey  fringe,  centered 

r1+r 

around  A  =  0,  agrees  closely  with  1/2 T  (-  628.9  Hz).  For  a  rectangular  shape  of  the 
Raman/CPT  pulse,  the  frequency  spread  of  the  FDRI  envelope  is  described  by  sine2 (At 7/2), 
which  corresponds  to  the  Fourier  transform  of  the  CPT  pulse.  Here,  we  have  used  the  function 
sine,  which  is  defined  as  sinc(x)=sin(x)/x. 


(a)  (b) 

Fig.  4.  Numerical  solutions  showing  (a)  FDRI  obtained  by  using  the  pulse  parameters  as 
described  in  the  text,  (b)  QPLD-RI  corresponding  to  different  values  of  A,  as  a  function  of  the 
query  pulse  duration. 

Figure  4b  shows  numerically  computed  QPLD-RI  fringes  for  different  values  of  A.  In  this 
case,  an  attenuated  (average  Rabi  frequency  27r.(T / 50) )  query  pulse  was  used,  with  z1  and  T 

chosen  to  be  212  ps  and  508  ps,  respectively.  Here,  the  horizontal  axis  is  the  duration  of  the 
query  pulse,  extending  to  a  value  exceeding  1  msec.  Thus,  these  oscillations  are  in  the  time 
domain.  The  time  period  of  these  oscillations  agree  with  1/A  ,  as  expected.  For  a  given  value 
of  A,  the  amplitude  of  the  fringe  now  depends  on  the  duration  of  the  query  pulse,  in  an 
oscillatory  manner.  Increasing  the  Rabi  frequency  (or  Raman  pumping  rate  a)  of  the  query 
pulse  causes  a  faster  decay  of  the  fringe  amplitude.  Negative  values  of  A  generate  the  same 
fringes  as  the  positive  values,  but  with  a  tz  phase  shift.  Thus,  for  a  give  value  of  large  x2,  the 
fringe  pattern  as  a  function  of  A  will  be  asymmetric,  especially  when  p^O.  In  addition,  the 
width  of  the  FDRI  fringes  will  depend  on  the  length  of  the  query  pulse.  The  asymmetry  and 
the  variation  in  the  width  can  both  serve  as  a  potential  source  of  error  for  atomic  clocks,  as 
discussed  in  greater  detail  later.  Next,  we  present  experimental  results  showing  these 
interference  effects  generated  in  a  buffer-gas  filled  rubidium  vapor  cell. 

3.  Experimental  setup 

Figure  5  shows  the  experimental  setup  used  to  study  QPLD-RI  and  FDRI  in  85Rb  atomic 
vapor.  A  10  cm  long  vapor  cell  containing  a  natural  mixture  of  Rb  isotopes  and  Ne  buffer  gas 
at  10  torr  pressure  is  used.  An  external-cavity  diode  laser  (ECDL,  Vortex  Series,  New  Focus) 
operating  at  795  nm  is  used  for  optical  excitation  of  the  D1  line  transitions  in  85Rb  atoms,  and 
is  locked  to  the  (F=3-F'=3)  transition  using  saturated  absorption  spectroscopy  (SAS). 

Part  of  the  ECDL  light  is  amplified  by  injection  locking  a  slave  laser  with  1-2  mW  seed 
power.  This  injection  locked  laser  (ILL)  uses  a  150  mW  diode  that  runs  near  800  nm  when 
not  injection-locked.  The  ILL  wavelength  is  monitored  using  an  optical  spectrum  analyzer  to 
ensure  single-wavelength  operation  with  >40  dB  sidemode  and  pedestal  suppression.  This 
amplified  light  passes  through  an  acousto-optic  modulator  (AOM)  centered  at  80  MHz  that 
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can  be  switched  on  with  <  1  (is  rise  time.  The  frequency  of  this  switched  light,  coh  is  also  at 
the  F=3  -  F'=3  transition  because  it  is  compensated  by  an  AOM  in  the  SAS  setup,  which  is  not 
shown  in  the  figure. 


Fig.  5.  Experimental  setup  for  Ramsey  interference  consisting  of  an  injection-locked  diode 
laser,  AOM  for  generating  Raman  beams  and  a  buffer  gas  filled  Rb  cell.  FCL:  fiber  collimating 
lens,  FCP:  fiber  coupler,  NDF:  neutral  density  filter,  OF:  optical  fiber,  P:  prism  (out-of-plane 
as  discussed  in  the  text),  PP:  prism  polarizer,  PBS:  polarizing  beam  splitter,  PD:  photodiode, 
and  HWP:  half-wave  plate. 

The  frequency- shifted  Raman  beam  at  co2  is  generated  from  co1  using  a  Bragg  acousto 
optic  modulator  (AOM;  f  =  1.5  GHz,  BW  =  ±100  MHz,  Brimrose  Corporation)  in  a  double¬ 
pass  configuration  as  shown  in  Fig.  5  for  operation  near  the  ground  state  hyperfine  frequency 
cohf  (=3.0357  GHz).  A  narrowband  microwave  voltage  controlled  oscillator  (VCXO)  driven 
near  1517.8  MHz  allows  us  to  vary  the  Raman  detuning  A  over  roughly  0.5  MHz  (100 
kHz/V).  The  input  light  for  the  AOM  is  fiber-delivered  and  focused  to  a  spot  size  of  60 
microns.  Since  the  AOM  diffraction  efficiency  is  polarization  dependent,  use  of  a  polarizing 
beam  splitter  at  the  input,  in  conjunction  with  a  polarization  rotation  after  the  first  pass,  for 
separating  the  double-pass  beam  from  the  input,  is  not  possible.  Instead,  The  shifted  beam 
after  the  first  pass  was  spatially  separated  from  the  input  beam  by  displacing  the  beam  out-of¬ 
plane  by  a  microprism  and  picked  off  by  an  edge  mirror.  The  out-of-plane  shift  causes  no 
noticeable  reduction  in  diffraction  efficiency,  and  the  total  double-pass  efficiency  of  the  AOM 
is  -5%.  The  undeflected  light  (&>;)  and  the  deflected  light  (co2)  are  each  coupled  into  separate 
polarization-maintaining  (PM)  fibers. 

These  two  beams  are  collimated  to  1/e  diameters  of  4mm  using  fiber  aspheres  (Thorlabs 
C260TM-B,  f=  15. 36mm)  with  near-Gaussian  profiles  and  orthogonal,  linear  polarizations  for 
Raman  excitation.  They  are  combined  using  a  polarizing  beam  splitter  for  copropagating 
excitation.  Both  beams  can  be  pulsed  simultaneously  using  the  AOM  in  the  ILL  setup,  and  A 
is  changed  by  sweeping  the  VCXO. 

The  cell  is  magnetically  shielded  and  heated  to  a  steady-state  temperature  near  50°  C  for 
producing  a  moderate  atomic  density  (1012/cm3).  Helmholtz  coils  inside  the  shielding  volume 
produced  a  small  homogeneous  axial  field  with  flux  density  <  100  mG.  This  causes  an 
antisymmetric  shift  of  the  Zeeman  sublevels  and  allows  us  to  separate  the  Raman  transitions 
between  |  F=2,  m,  =0)  and  |  F=3,  m,  =  0)  (the  so-called  clock  transition)  from  the  other, 

magnetically- sensitive  Raman  transitions.  The  clock  transition  is  insensitive  to  magnetic  field 
to  a  first-order,  but  one  can  also  choose  magnetically- sensitive  Raman  transitions  if  interested 
in  doing  sensitive  magnetometry. 
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A  high  extinction  ratio  Glan-Thompson  polarizer  (40  dB)  is  used  after  the  vapor  cell  to 
remove  one  of  the  Raman  beams  from  reaching  the  detector.  Ramsey  interference  in  the  pulse 
amplitude  is  observed  by  a  gated  boxcar  integrator  [14]  which  generates  an  integrated  signal 
corresponding  to  the  amplitude  of  the  query  pulse.  Note  that  Ramsey  interference  can  also  be 
observed  without  the  prism  polarizer  since  the  effect  is  manifested  in  both  Raman  beams.  The 
prism  polarizer  allows  us  to  observe  the  effect  at  one  frequency  (cu7  or  co2)  without  being 
encumbered  by  the  beating  with  the  other. 

4.  Experimental  results  and  discussions 

Two  different  pulsed  Raman  excitation  schemes  are  used  in  the  experimental  setup  shown  in 
Fig.  5  to  generate  QPLD-RI  and  FDRI  in  85Rb  vapor.  In  both  cases,  relatively  long  CPT 
pulses  (zj  =  200  ps)  are  used  for  the  initial  Raman  excitation.  FDRI  is  observed  by  sending  a 
repeated  sequence  of  Raman  pulse  pairs  (generated  by  the  AOM  in  the  ILL  setup)  consisting 
of  a  short  query  pulse  (r2  =  400  ns),  and  by  scanning  the  frequency  of  co2  to  change  the  two- 
photon  resonance  condition.  The  pulse  sequence  is  synchronized  with  the  voltage  ramp  that 
controls  the  VCXO. 


0.4 


Fig.  6.  (a)  Experimentally  observed  FDRI  using  Raman  excitation  in  the  (0-0)  ground-state 
‘clock  transition’  in  85Rb  vapor,  (b)  second-order  magnetic  field,  B  dependence  of  FDRI. 


Figure  6a  shows  experimentally  observed  FDRI  using  T  =  800  ps.  The  result  qualitatively 
matches  our  numerical  simulation  shown  earlier.  The  line  width  Av  is  inversely  related  to  the 
free  evolution  time  T  between  the  Raman  pulse  pair  and  is  measured  to  be  close  to  625  Hz.  A 
total  optical  power  of  about  300  pW  is  used  in  the  Raman  beams  and  a  power  ratio  of  2:1 
between  the  Raman  beams  is  used  in  this  case,  in  order  to  ensure  that  the  Rabi  frequencies  for 
the  two  legs  of  the  A-transition  are  nearly  equal,  taking  into  account  the  different  oscillator 
strengths  along  the  two  legs.  In  our  experiments,  T  of  up  to  1.5  ms  has  been  used  to  generate  a 
fringe  width  as  narrow  as  330  Hz,  with  a  signal-to-noise  ratio  of  about  10.  The  maximum 
value  of  T  is  limited  by  the  ground- state  decoherence  time  in  the  medium. 

The  central  frequency  in  the  above  experiment  is  defined  by  the  |  F=2,  mF=0)  and 

|F=3,mF=0)  Raman  transition  in  85Rb  atoms.  It  has  a  quadratic  (second-order)  dependence 

on  the  magnetic  field.  The  sensitivity  of  this  frequency  to  the  external  magnetic  field,  B  is 
verified  by  studying  the  magnetic  field  dependence  of  FDRI.  This  is  done  by  changing  the 
strength  of  axial  magnetic  field,  B,  linearly  using  the  Helmholtz  coil,  and  keeping  the 
frequency  difference  between  the  Raman  beams  matched  to  the  condition  A  -  0  for  B=0. 
Figure  6b  shows  FDRI  fringes  for  T  =  400  ps,  while  scanning  B  around  a  dc  (or  bias) 
magnetic  field  (approximately  2.5  G)  so  that  it  does  not  change  its  direction.  The  frequency  of 
the  interference  pattern  is  chirped  due  to  the  quadratic  dependence  on  B.  This  can  be 
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described  [23]  as  /(#)  -  (l-e  ^  )-[l  +  cos (a0_0  B2T  ,  where  77 


r  % 

1  +  3^  ’ 


«o-o  (  =  1-293 


kHz/G2)  is  second-order  magnetic  field  coefficient  associated  with  the  clock  transition,  while 
assuming  that  ij  »  r2. 

QPLD-RI  was  observed  in  the  experiment  by  replacing  the  short  detection  pulse  used  in 
FDRI  by  a  long,  attenuated  query  pulse  of  duration  up  to  z2  -  4  ms  in  the  Raman  pulse 
sequence.  In  this  case,  Ramsey  interference  fringes  were  observed  directly  as  time-domain 
oscillations  in  the  amplitude  of  the  query  pulse  without  the  use  of  boxcar  integrator.  The 
optical  power  in  the  query  pulse  was  attenuated  by  nearly  twenty  times  so  that  it  does  not 
cause  any  significant  decay  of  coherence  created  by  the  first  Raman  pulse.  Figure  7a  shows 
QPLD-RI  observed  in  the  pulse  amplitude  due  to  Ramsey  interference,  as  explained  before. 
The  vertical  axis  for  each  plot  corresponding  to  different  values  of  A  is  displaced  in  order  to 
show  all  the  plots  together.  The  initial  part  of  the  plot  also  shows  the  first  Raman  pulse 
preceding  the  query  pulse.  For  this  experiment,  z1  and  T  were  chosen  to  be  400  ps  and  250  ps, 
respectively.  The  period  of  these  oscillations  closely  match  with  1/A  .  The  decay  of  the 


fringes  is  due  to  a  combination  of  two  effects:  saturation  of  the  system  into  the  dark  state,  and 
decoherence  due  to  collision.  For  our  experiment,  the  decoherence  is  much  faster  than  the  rate 
at  which  the  system  decays  to  a  dark  state.  Under  this  condition,  the  number  of  observed 

oscillations  in  QPLD-RI  can  be  described  as  /(r)  -  |^l  +  cos(A.  (r  +T)  )Jexp(-£/Q) ,  where 

tc  is  the  ground-state  decoherence  time.  Our  results  indicate  that  the  coherence  in  the  medium 
survives  approximately  for  3-4  ms.  This  illustrates  the  fact  that  the  coherence  survives  many 
collisions  as  the  atoms  slowly  diffuse  in  the  interaction  region  in  the  presence  of  the  buffer 
gas  [25]. 
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Fig.  7.  (a)  Experimentally  observed  QPLD-RI  for  different  values  of  A  as  shown  in  the  figure, 

(b)  magnetic  field  dependence  of  QPLD-RI  for  A  «  0.  The  horizontal  axis  shows  the  whole 
pulse  sequence. 

Transient  coherent  oscillation  has  been  reported  by  making  changes  to  the  Raman 
detuning  for  CPT  fields  [26].  We  would  like  to  point  out  that  the  effect  reported  in  ref.  26  is 
quite  different  from  QPLD-RI.  The  oscillations  reported  there  was  observed  while  changing 
the  Raman  detuning  as  a  function  of  time,  either  abruptly  or  continuously.  Furthermore,  there 
was  no  dark  time  interval.  In  contrast,  QPLD-RI  is  observed  while  the  Raman  detuning  is  kept 
constant,  and  there  is  a  dark  time  interval  between  the  initial  excitation  and  the  query  pulse. 
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We  have  also  studied  the  magnetic  field  dependence  of  QPLD-RI  by  applying  different 
axial  magnetic  fields.  Figure  7b  shows  QPLD-RI  obtained  by  using  5=  1.5  G,  while  keeping 
A  -  0.  The  pulse  parameters  chosen  in  this  case  are  r7  =  400  ps,  T  =  250  ps  and  r2  -  4  ms. 

These  oscillations  can  be  described  as  /(r)  -  ^l+cos(a0_0B2  +  )Jexp(-f  / tc) .  If  a 

different  magnetic  field  is  applied,  the  periodicity  of  the  fringes  changes  quadratically  (not 
shown)  with  the  field  strength.  The  effective  decoherence  time  measured  from  the  decay  of 
the  fringe  amplitude  is  shorter  (tc  -  2  ms)  compared  to  those  shown  in  Fig.  7a.  This  is 
possibly  due  to  the  inhomogenity  of  the  axial  magnetic  field  over  the  length  of  the  cell. 

Finally,  we  discuss  how  the  QPLD-RI  can  affect  the  accuracy  of  an  atomic  clock  under 
various  conditions.  Consider  first  the  ideal  case  where  the  average  optical  detuning  3  is  set  to 
zero.  Figure  8  shows  the  effect  of  choosing  different  query  pulse  duration,  x2,  on  the  FDRI 
fringes,  using  numerical  calculations.  The  Rabi  frequencies  for  Raman  beams  were  chosen  to 
be  =  Q2  =  2k  (r/5)  during  the  first  Raman  pulse,  and  f27  =  Q2  =  2n  (r/200)  during  the 
second  (or  query)  Raman.  We  used  a  repeated  pulse  sequence  in  the  calculation,  with  pulse 
parameters  r7  «212  ps  and  T  =  795  jus ,  same  as  shown  earlier.  Here,  we  show  only  the  central 
FDRI  fringe  (Fig.  8a),  for  three  different  query  pulse  durations:  10  ps,  100  ps,  and  300  ps. 


Fig.  8.  Numerical  results  showing  change  in  (a)  central  FDRI  fringe  width  and  (b)  slope  at  the 
fringe  center  as  a  function  of  query  pulse  duration,  x2. 

In  Fig.  8b,  we  have  plotted  the  central  slope  of  the  derivative  of  the  fringes,  for  a  range  of 
query  pulse  durations,  r2.  This  slope  is  a  measure  of  the  width  of  each  fringe.  As  can  be  seen, 
the  slope  increases  with  query  pulse  length,  corresponding  to  narrowing  of  fringes.  It  is  well 
known  that  the  shot-noise  limited  fractional  frequency  stability  (FFS)  of  a  clock,  stabilized  to 
the  center  of  one  such  fringe,  is  given  by  (j  /  /0)  /  (5  /  N) ,  where  y  is  the  fringe  width,  f0  is 

the  central  frequency  of  the  clock  (matching  the  frequency  difference  between  I3>  and  ll>), 
and  (5/AO  is  the  signal-to-noise  ratio  (SNR).  Thus,  for  a  given  SNR,  the  FFS,  which  is  the 
defining  parameter  for  a  clock  used  in  any  system,  will  vary  if  the  query  pulse  length  is 
varied.  This  dependence  is  currently  not  taken  into  account  in  considering  the  FFS  of  a  pulsed 
Raman  clock.  A  key  message  of  this  paper  is  that  one  must  pay  attention  to  the  length  of  the 
query  pulse  in  determining  the  FFS  of  such  a  clock. 

Consider  next  the  non-ideal  case  where  the  average  detuning,  3 ,  is  non-zero.  In  any 
Raman  clock,  this  condition  is  relevant  in  order  to  account  for  inevitable  fluctuations  in  the 
laser  frequency.  For  a  vapor  cell  base  Raman  clock,  this  condition  is  also  unavoidable  because 
of  the  presence  of  several  hyperfine  transitions  in  the  excited  manifold  that  can  all  contribute 
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to  the  clock  signal,  due  to  Doppler  broadening.  For  illustrative  purposes,  we  first  choose  8  =  - 
T  in  Fig.  9a.  Here,  we  show  the  FDRI  fringes  for  three  different  query  pulse  durations:  10  ps, 
210  ps,  and  300  ps.  Note  that  in  each  case  the  center  of  the  fringe  peak  is  shifted  from  the 
ideal  position  of  A  =  0.  This  shift  is  proportional  to  the  intensity  of  the  fields  used  in  the  first 
excitation  pulse,  and  is  thus  a  manifestation  of  the  light-shift  or  the  ac-Stark  shift  [23].  We  see 
that  as  the  query  pulse  length  increases,  the  amount  of  shift  changes,  and  the  fringes  become 
increasingly  asymmetric.  In  Fig.  9b,  we  have  plotted  the  shift  as  a  function  of  the  query  pulse 
duration,  for  8  =  -1.5 T  .  In  this  case,  the  variation  is  seen  to  be  approximately  linear;  in 
general,  the  exact  functional  dependence  of  the  shift  on  the  query  pulse  is  expected  to  be 
different  for  different  choice  of  parameters.  Currently,  the  laser  line  width,  the  Doppler  width, 
the  relative  frequency  separation  between  the  hyperfine  levels  in  the  excited  state  manifold, 
and  the  dependence  of  light  shift  on  8  together  are  used  to  estimate  the  light- shift  induced 
limit  of  the  clock  accuracy  for  a  vapor  cell  based  Raman  clock.  From  the  result  illustrated  in 
Fig.  9,  we  conclude  that  the  query  pulse  duration  must  also  be  taken  into  account  in 
determining  this  limit  accurately. 
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Fig.  9.  Numerical  results  showing  (a)  light  shift  in  central  fringe  for  5  =  -T,  and  (b)  variation  in 
light  shift  as  a  function  of  x2  corresponding  to  5  =  -1.5r. 

The  results  presented  in  this  paper  are  only  qualitative.  For  a  quantitative  determination  of 
the  effect  of  query  pulse  dependence  for  a  vapor  cell  based  Raman  clock,  one  must  take  into 
account  the  Doppler  shift  for  atoms  moving  in  three  dimensions,  the  hyperfine  levels  in  the 
excited  state  manifold,  the  different  strengths  of  the  various  transitions,  and  optical  pumping 
into  and  out  of  various  Zeeman  sublevels.  We  are  in  the  process  of  carrying  out  such  an 
analysis.  Until  the  comprehensive  theoretical  modeling  is  completed,  we  cannot  expect  to 
predict  precisely  the  effects  of  QPLD-RI  on  the  clock  frequency  fluctuations.  Furthermore,  in 
our  clock  setup,  excess  noise  induced  fluctuations  are  too  large  to  observe  these  effects 
currently.  Efforts  are  underway  to  improve  our  apparatus  to  reach  shot-noise  limited 
performance,  so  that  these  effects  can  be  measured.  We  expect  to  report  on  the  comprehensive 
theoretical  model  as  well  as  experimental  studies  of  these  effects  in  the  near  future. 

5.  Conclusion 

In  summary,  we  have  shown  a  new  type  of  Ramsey  interference  effect,  QPLD-RI,  which  is 
observed  when  a  long,  attenuated  query  pulse  is  used  during  pulsed  Raman  excitation  in 
rubidium  vapor.  This  effect  can  produce  pulse  length  dependent  broadening  as  well  as 
asymmetry  in  FDRI  fringes  used  for  a  vapor  cell  Raman  clock,  causing  a  potential  source  of 
error.  Specifically,  we  have  shown  that  broadening  affects  the  fractional  frequency  stability  of 
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a  clock,  and  asymmetry  changes  the  magnitude  of  the  light-shift  induced  error  due  to  optical 
detuning,  which  is  unavoidable  in  a  vapor  cell  based  Raman  clock  because  of  the  presence  of 
several  hyperfine  transitions  in  the  excited  state  manifold,  which  all  contribute  to  the  clock 
signal  because  of  Doppler  broadening.  Thus,  in  choosing  the  optimal  parameters  for  such  a 
clock,  the  effect  of  QPLD-RI  has  to  be  taken  into  account.  The  origin  of  FDRI  and  QPLD-RI 
has  been  explained  qualitatively  using  a  Bloch  vector  model  as  well  as  numerical  solutions  to 
time-dependent  density  matrix  equations. 
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